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1 Introduction 

Recent neutron scattering experiments by Aeppli et aM have studied the two- 
dimensional incommensurate spin correlations in the normal state of the high 
temperature superconductor La a; Sr 1 _, I .Cu04 at x « 0.15 in some detail. In 
particular, they have determined the scattering cross section over a signifi- 
cant portion of the wavevector, frequency and temperature space. One of 
their striking observations has been that, while the dynamic structure factor 
of the electronic spins has quite an intricate functional form over this three- 
dimensional parameter space, the results become quite simple and explicable 
when interpreted in terms of the non-zero temperature dynamical properties 
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of a system in the vicinity of a quantum critical point au. The data are consis- 
tent with such an interpretation over a decade in temperature and in over two 
decades of the peak scattering amplitude, and indicate that there is a nearby 
quantum critical point with dynamic critical exponent The nature of the 

measured spin correlations indicates that this quantum critical point is the po- 
sition of a quantum phase transition to an insulating, charge- and spin-ordered 
ground state. Varying the doping concentration, x, alone is not expected to 
be sufficient to access this ordered state, and authors El asserted that a second 
tuning parameter ('j/') is necessary; e.g. it is known that replacing some of the 
Sr by Nd allows one to move along the y axiscl. 

In this paper, we will review recent theoretical work studying non-zero 
temperature (T) dynamical response functions closely related to those mea- 
sured in the neutron scattering experiment of Aeppli et al. Ill We will examine 
a simple class of models in spatial dimension d = 1,2, which exhibit quantum 
ordering transitions with z = 1. All of the models we shall study are closely 
related to the following quantum field theory of a real, n-component, scalar 
field 4> a (a = 1 . . . n): 



Z Q = J V<t> a (x,T)exp (- J d d x 



l/T 

drCc 



Cq -\ 



^(^ Q ) 2 + (V^ Q ) 2 + (r c + r)^ +^(c (11 



Here we are using units with h = fcs = 1, x is the d-dimensional spatial co- 
ordinate, r is imaginary time, c is a velocity, and r c , r and u are coupling 
constants. The co-efficient of the 4>a term (the 'mass' term) has been written 
as r + r c for convenience. The quartic non-linearity u makes Zq as interacting 
quantum field theory, and it ultimately responsible for the dynamic relaxation 
phenomena we shall describe. For appropriate values of d and n, the model Zq 
exhibits a T = quantum transition between an ordered phase with (<p a ) =/= 0, 
and a quantum paramagnet with complete O(n) symmetry; the value of r c 
is chosen so that this transition occurs at r = 0. In physical applications 
of Zq, the case n = 1 describes the Ising model in a transverse field, the 
case n = 2 case describes a superfluid-insulator transition (with cf>i + ifo, the 
complex superfluid order parameter), and n — 3 describes spin fluctuations 
in a collinear quantum antiferromagnet (with (f> a the antiferromagnetic order 
parameter) . 

We will be interested in the real-time dynamic properties of Zq in the 
high T limit of the continuum, quantum field theory (in some cases, this is also 
the 'quantum critical' region a). More specifically, in the vicinity of a quantum 
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critical point, the theory Zq can be characterized by two distinct energy scales. 
The first, is a low energy scale, characterizing the deviation from the system 
from the quantum critical point at r = 0: this varies as 6|r| 2ly , where v is 
the correlation length exponent, and b is a non-universal, cutoff-dependent 
constant needed to make the expression have physical units of energy; this low 
energy scale is the central quantity characterizing the continuum quantum field 
theory. The second, is a high energy scale, and is of order cA, where A is a 
high-momentum cutoff needed at intermediate stages to define the continuum 
limit of Zq; this energy scale plays no role in the continuum quantum field 
theory. We will be assuming here that the temperature is in between these two 
energy scales i.e. 

b\r\ zlJ « T « cA. (2) 

We will characterize the response of the system by the dynamic suscepti- 
bility, x(k,u n ), defined by 

X (k,u; n ) = - [ 1/T dr f d d x V (<(> a (x, r)4> a (0, 0)} e -^-^r) ; (3) 

where k is the wavevector, u> n the imaginary frequency; throughout we will use 
the symbol uo n to refer to imaginary frequencies, while the use of u will imply 
the expression has been analytically continued to real frequencies. Further, 
the static susceptibility, x(fc), is defined by 

X(k) = x (k,io n = 0). (4) 

We shall be especially interested here in a relaxation function 1Z(k,u), which 
we define as 

1 2 

K(k,uj) = ——-lmx(k,w). (5) 
X(k) v 

This is an even function of u> with the dimensions of time, and the Kramers- 
Kronig relation implies that 

/OO 7 
— ft(fc, W ) = l. (6) 
-oo 27T 

After a Fourier transform to real time, 1Z(k,t) (with lZ(k,t = 0) = f) de- 
scribes the time-dependent relaxation of spin correlations at wavevector k due 
to quantum and thermal fluctuations. 

(In a regime where the predominant fluctuations have an energy much 
smaller than T, the low frequency dynamics becomes effectively classical, and 
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then the fluctuation-dissipation theorem implies that 



where S(k,uj) is the dynamic structure factor, and S(k) is the equal-time 
structure factor. However, this relationship is not generally true for a quantum 
system, and for clarity, we will always use (||) as our defining relation.) 

The following sections will describe the behavior of TZ(k,u>) in the high T 
limit of the continuum quantum theory Zq for various physically interesting 
cases in low dimensions: we will consider the case n = 1, d = 1 in Section 
the case n = 3, d = 1 in Section ||, and all values of n in d — 2 in Section ||. 

2 Ising chain in a transverse field 

The Ising chain in a transverse field has the Hamiltonian 

where erf ' 2 are Pauli matrices on the sites, i, of a d = 1 chain with lattice 
spacing a. This model has a quantum critical point at g = 1, whose vicinity 
is believed to be described by the d = 1, n = 1 case of Zq. Under this 
mapping, the field-theoretic correlators of (j) a under Zq are equivalent to the 
long-distance, long-time correlators of the order parameter a z under the lattice 
model Hj. 

The continuum high T dynamical response of (^) can be computed exactly 
by a special trick relying on the conformal invariance of the critical theory at 
g = lu. The result for x(k,u) is 

Zc G/(0)r(7/8) r 
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X(A,W) ^ ^ p^g^ /IK ,. , „L\ /1C .. „;„\ ■ ( 9 ) 



Here Z = J" 1 / 4 , c = 2Ja, and G/(0) = 0.858715 .... We show a plot of Im* in 
Fig 1. For uj,ck^> T there is a well-defined 'reactive' peak in Im% at uj w c/c 
reflecting the excitations of the T = quantum critical point, which have an 
energy threshhold at 10 = ck. However the low frequency dynamics is quite 
different, and for uj,ck «T we cross-over to the quantum relaxational regimeul. 
This is made clear by an examination of Imx(fc, uj)/uj as a function oiuj/T and 
cfc/T, which is shown in Fig 2. Now the reactive peaks at uj ~ ck are just about 
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Figure 1: High temperature dynamic susceptibility, T 7 / 4 Imx(fc, oj)/Z in (|9p of the quantum 
Ising chain [d = 1, n = 1) as a function of lo/T and ck/T. 



invisible, and the spectral density is dominated by a large relaxational peak 
at zero frequency. We can understand the structure of Fig 2 by expanding the 
inverse of (|^) in powers of k and u>; this expansion has the form 

X(k,u) = , (10) 

1 - i(w/ui) + k^e - (W/W2) 3 

where x(0) ~ T -7 / 4 , and lj\.2 and £ are parameters characterizing the expan- 
sion. For k not too large, the w dependence in (p!o|) is simply the response of 
a strongly damped harmonic oscillator: this is the reason we have identified 
the low frequency dynamics as "relaxational" . The function in ( [To| ) provides 
an excellent description of the spectral response in Fig 2. We determined the 
best fit values of the parameters lo\^ and £ by minimizing the mean square 
difference between the values of Imx(k,u>)/u> given by @ and @ over the 
range < u> < 2T and < ck < 2T and obtained 

ui = 0.396 T 
u 2 = 0.795 T 
1= 1.280 c/T. (11) 
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Figure 2: Plot of the spectral density of the quantum Ising chain (d = 1, n = 1) 
T 11 f A lvcix(k,uj)/ujZ as a function of u>/T and ck/T. Note that this is simply the quan- 
tity in Fig 1 divided by u). The reactive peaks at u) ~ ck in Fig 1 are essentially invisible, 
and the plot is dominated by a large relaxational peak at zero wavevector and frequency. 



The quality of the fit is shown in Figs 3 and 4. In Fig 3 we compare the 
predictions of and ( fiol) for lmx(k , lo) / 'lo at lo = as afunction oick/T. The 
form ([To]) predicts a Lorentzian-squared response function and this is seen to 
provide a better fit than a Lorentzian — a similar Lorentzian-squared response 
was used in analyzing the data in Ref . 1 . In Fig 4 we plot the predictions of (||) 
and © for TZ(k, lo) at ck/T = 0, 1.5 as a function of uj/T. For k = (lo = 0) 
there is a large overdamped peak at u> — (k = 0), but a weak reactive peak 
at to ~ ck does make an appearance at larger wavevectors or frequencies. 

For an alternative, and more precise, characterization of the relaxational 
dynamics we can introduce the relaxation rate Tr defined by 

r;' , » ; (i2) 



we have chosen this definition because for the suggestive functional form (Q) , 
Tr = loi, the frequency characterizing the damping. However, using (||) we 
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Figure 3: Comparison of the predictions of (Kj) (dots) and (gxj) (solid line) for lmx{k,u>) /u> 
of the quantum Isins chain (d = 1, n = 1) atjji = as a function of ck/T. The best 
fit parameters in (hll) were used. The function (jlOh yields the square of a Lorentzian as a 
function of k; a best fit by just a Lorentzian is also shown (dashed line), and is much poorer. 



determine: 

/ vr \ k B T 

r - = ( 2tan i6)i- 

7 rr\ 

w 0.397825-1—, (13) 
n 

where we have inserted physical units to emphasize the universality of the 
result. 

An important property of the present high T dynamical results is that at 
the scale of the characteristic rate T^, the dynamics of the system involves in- 
trinsic quantum effects which cannot be neglected. Description by an effective 
classical model would require that Tr <C ksT/U, which is thus not satisfied 
here. 

To obtain an intuitive physical picture of the above quantum relaxational 
dynamics, it is useful to consider approaching the high T limit by gradually 
rising the temperature (while keeping all other couplings fixed) from the low 
T limit where T < b\r\ zy . 
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First, consider the low T limit on the magnetically ordered side. Here the 
excitations above the ground state are 'domain walls' which separate regions 
in which the Ising order parameter has opposite signs. These domain walls can 
move easily without significant change of energy, and their low energy motion 
leads to a large relaxational peak in lZ(k, uj) at uj = 0, k — 00. At very low T, 
the domain walls are very dilute, and their spacing is much larger than their 
thermal de Broglie wavelengths — consequently their motion can be described in 
a classical model. However, as T is raised into the high T regime, their spacing 
becomes of order their de Broglie wavelength, and the relaxation rate of their 
collisional dynamics becomes of order T: this is leads to the relaxational peak 
in Fig 4. 

Second, we can begin by considering the low T limit on the quantum 
paramagnetic side. Now the excitations are local 'flipped spins' which require 
a finite energy, A, to create them. So there is a sharp peak in lZ(k,uj) at 
w = A, which is broadened by collisions with the dilute, classical gas of pre- 
existing quasiparticles. In the language of the field <p a , this finite frequency 
peak arises from amplitude fluctuations in cf> about a local minimum in its 
effective potential. As T is raised, the quasiparticle gas becomes dense with 
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the mean-particle spacing becoming of order their de Broglie wavelength, the 
quasiparticle line- width becomes of order T, and the peak in lZ(k, uj) eventually 
moves to lo = 0. 



3 Gapped, Heisenberg antiferromagnetic chains 

In this section we will consider Zq for the case d = 1, n = 3, which describes the 
low energy fluctuations of certain one-dimensional quantum antiferromagnets: 
chains of integer spin S, or p-leg ladders with p even. The ground state of these 
systems is always a quantum paramagnet with an energy gap, A, and there 
is no magnetically ordered state. So strictly speaking, there is no quantum 
critical point, and it may appear that the general arguments of Section |l| do 
not apply. However, it is still possible to define a universal, continuum high T 
limit of the quantum theory. We have to replace the condition ^ by 

A < T < cA ~ J, (14) 

where J is a typical exchange constant of the antiferromagnet, but, unlike (0), 
the energy gap, A, is always non-zero. So we have to pick a system with A 
much smaller than any microscopic energy scale. As A becomes exponentially 
small as S (or p) is increased, this is quite easily achieved. In a formal sense, we 
can consider the following an analysis of the high T dynamics of the 'quantum' 
critical point reached in the limit S — ► oo, when A becomes vanishingly small. 

The high T limit introduced above has been analyzed in some detail in 
two recent papers uu. It was argued that there is an effective classical non- 
linear wave problem that describes the long time relaxation dynamics. The 
degrees of freedom of this classical model are a 3-component unit length field 
n(x,t), n 2 = 1, which describes the local orientation of the antiferromagnetic 
order parameter, and its canonically conjugate angular momentum, L(x,£), 
which measures the ferromagnetic component of the local spin density. The 
equal time correlations of L and n are described by the following continuum 
classical partition function 

Hie 



T 



Zic = J Vn(x)VL(x)S(n 2 - 1)6 (L ■ n) exp 

(15) 

, ax J XuX 

The parameters in this partition function are determined universally jp. terms 
of A, T, and c, and exact expressions are known in the limit Q M. The 



, dn \ 2 1 , 
T£ — + L 2 
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correlation length, £, is given by 



^2^T ln ( A J' < 16 > 

where 7 is Euler's constant. The quantity Xu± is the susceptibility to a uniform 
magnetic field (which couples to the ferromagnetic moment) in a direction 
orthogonal to the local antiferromagnetic order; it is related to the rotationally 
averaged uniform susceptibility, \u, by 

2 

Xu = ^Xu±, (17) 

and Xu is given by 

1 , /32^H 2 +^T\ 
** = 3^H A J' ^ 

As we will see below, with these parameters in hand, the characteristic 
excitation of the classical model ( |l5| ) has energy s ~ T j ln(T/A). For T A, 
this is parametrically smaller than T. So the occupation number of the wave 
modes with energy e will be much larger than unity, and the quantum Bose 
function will take the classical equipartition value T/e. This is the argument 
which justifies use of a classical model in this high T limit. 

All equal time correlations of the model ([if]) can be computed exactly. As 
this is a model to which the classical fluctuation-dissipation theorem applies, 
the equal time, two-point n correlator is directly related to the static sus- 
ceptibility; the underlying quantum fluctuations however do induce an overall 
wavefunction renormalization factor Bia. The two point n correlator decays 
exponentially on the scale £, and by its Fourier transform to momentum space 
we obtain 



T x {k) = A 



In 



2{/3 P.) 



(l + fc 2 £ 2 )' 

Here A is a non- universal amplitude which determines the scale of the field <f> a , 
and the multiplicative logarithmic factor comes from the underlying quantum 
fluctuations; the remaining is just the Fourier transform of e~l x l/£/3, the 1/3 
coming from the 1/n in (|3|) . 

Let us now turn to the unequal time correlations. To obtain these, 
have to supplement Z\c with equations of motion, which have been arguedt 
to be the Hamilton-Jacobi equations associated with the Poisson brackets of 
n and its canonically conjugate angular momentum L: 

{L a (x), Lf}(x')} PB = e a /3-,L 7 (x)8(x - x') 



10 



{L a (x), np(x')} pB = e a i3 1 n 1 {x)5(x - x') 
{n a (x),n/3(x')} PB = 0. (20) 



From this, and (|15|) , we obtain directly the equations of motion for the quasi- 
classical waves 

— = {n,Hic}pB 
-L x n 



— = {L,Hic}pb 

dx 



(TOnx^. (21) 



To compute the needed unequal time correlation functions, pick a set of initial 
conditions for n(x), L(a;) from the ensemble (^). Evolve these deterministi- 
cally in time using the equations of motion (|^ ) . The value of the correlator is 
then the product of the appropriate time-dependent fields, averaged over the 
set of all initial conditions. We also note here that simple analysis of the dif- 
ferential equations (|2l|) shows that small disturbances about a nearly ordered 
n configuration travel with a characteristic velocity c(T) given by 

c(T) = {T^T)/ X ux{T)) l '\ (22) 

which is a basic relationship between thermodynamic quantities and the ve- 
locity c(T). Notice from ( [l6| ) and (JlJ) that to leading logarithms c(T) sa c, 
but this result is not satisfied by the subleading terms. The characteristic ex- 
citation will have energy e ~ c/£, and this leads to our estimate for e made 
earlier, when we justified the validity of a classical model. 

The classical dynamics problem defined by ( |l5| ) and (21) obeys that the 



crucial property of being free of all ultraviolet divergences. Consequently, we 
may determine its characteristic length and time scales by simple engineering 
dimensional analysis, as no short distance cutoff scale is going to transform 
into an anomalous dimension. Indeed, a straightforward analysis shows that 
this classical problem is free of dimensionless parameters, and is a unique, 
parameter-free theory. This is seen by defining 

_ x 

* = ? 
__ t 
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where the characteristic time, t v , is given by 



(24) 

our notation suggests that t v (like T^ 1 earlier) is a phase coherence time be- 
yond which rclaxational dynamics of damped spin waves takes over. Inserting 
(^Hi]) into ([is]) and ( pl| ) , we find that all parameters disappear and the par- 
tition function and equations of motion acquire a unique, dimensionless form, 
given by setting T — £ — Xu± = 1 in them. 

The above transformations allow us to easily obtain a scaling form for the 
relaxation function 1Z: 



-R{k,u)^T v ^ n {k^uoT v ), (25) 

where is a universal scaling function, normalized as in (^). Further infor- 
mation on the structure of ^>n was obtainedu by a combination of analytic and 
numerical methods. At sufficiently large we expect a pair of broadened, 
reactive, 'spin-wave' peaks at to « c(T)k (with c(T) given in (|2^)), which are 
similar to those found in the high T limit of the quantum Ising chain in Fig 4. 
For the opposite limit of small we present numerical results for ^h{Q,uj) 
in Fig 5. There is a sharp relaxational peak at uj = 0, which is again similar to 
that found in the high T limit of the quantum Ising chain in Fig 4. However, 
there is now a well-defined shoulder at uj w 0.7 which was. pot found in the 
Ising case. This shoulder is a remnant of the large n result at3 which predicts 
a delta function at w ~ T/ln(T/A). So N — 3 is large enough for this finite 
frequency oscillation to survive in the high T limit. 

There is alternative, helpful way to view this oscillation frequency. The 
underlying degree of freedom in our dynamical field theory has a fixed ampli- 
tude, with |n| = 1. However, correlations of n decay exponentially on a length 
scale £ — so if we imagine coarse-graining out to £, it is reasonable to expect 
significant amplitude fluctuations in the coarse-grained field. It is now useful 
to visualize an effective field <p a with no length constraint, which is just the 
field we introduced in Section fil. On a length scale of order £, we expect the 
effective potential controlling fluctuations of (f> a to have minimum at a non- 
zero value of \4>a\, but to also allow fluctuations in \<fi a \ about this minimum. 
The finite frequency in Fig 5 is due to the harmonic oscillations of ip a about 
this potential minimum, while the dominant peak at uj = is due to angular 
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fluctuations along the zero energy contour in the effective potential. This is 
interpretation is also consistent with the large n limit, in which we freely in- 
tegrate over all components of n, and so angular and amplitude fluctuations 
are not distinguished. The above argument could also have been applied to 
the quantum Ising chain (in this case, angular fluctuations are replaced by 
low-energy domain wall motion), but the absence of such a reactive, finite fre- 
quency peak at k = in Fig 4 indicates that n = 1 is too far from n = oo for 
any remnant of this large n physics to survive. 

4 Order parameter dynamics in d = 2 

Determination of the long time dynamics in d — 2 is a strongly coupled problem 
and remains unsolved. As in the analysis in Section |^, we expect a phase co- 
herence time, Tfi, and an inverse relaxation rate, T^ 1 , of order T, and quantum 
and thermal fluctuations to play an equally._important role. 

However, it has been argued recently til, that within the context of an 
expansion in 

e = 3 - d, (26) 
it is possible to separate quantum and thermal effects into two distinct stages 
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of the calculation, and to derive an effective classical wave model for the long 
time dynamics. As in Section]^, we first derive an effective action for the static 
susceptibility, and then supplement it with equations of motion to obtain the 
unequal time correlations. We define the zero Matsubara frequency component 
of (j) a by 



and derive an effective action for <I> Q by integrating out the non-zero Matsubara 
frequency components of <p a - To leading non-trivial order in an expansion in 
e, this leads to the following effective action 



As in (|15|), along with the functional integral over <& a (x), we have included 
an integral over a conjugate momentum field Tl a (x) which will be important 
for our subsequent treatment of the dynamics; for now it easy to see that 
the Gaussian integral over II Q can be performed exactly, and it leaves the 
correlations of <I> Q under Z 2( j unchanged. The coupling constants in ( p8[ ) are 
universally related to the underlying field theory Zq controlling the quantum 
critical point. Before specifying these, it is crucial to understand the nature of 
the ultraviolet divergences in Z 2 c considered as a_classical field theory in its 
own right. From standard field-theoretic analyses Ej it is known that Z 2 c has 
only one ultraviolet divergence, coming from a single one-loop tadpole graph in 
the self energy: consequently, by trading the bare 'mass' R for a renormalized 
mass R defined by 



we can remove all cutoff dependencies in the correlators of Z 2 c (there are 
some additional divergences, associated with composite operators, which ap- 
pear when two or more field operators approach each other in space: we will 
not be concerned with these here). All observables of Z 2 c are then universal 
functions of the couplings R and U. Moreover, it is precisely these couplings 
that are universally computed from the underlying quantum field theory Zq. 
Actually, instead to dealing with R and U as the two independent parameters 




(27) 




(28) 




(29) 
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controlling correlators of Z2C, it is convenient to replace U by the dimension- 
less parameter Q denned by 



TU 



R (i-d)/2-' 



(30) 



which is analogous to the Ginzburg parameter-, ia the high T limit (g), these 
parameters were shown to have the following universal values to leading 
order in e 



R = e 



+ 2\ 2n 2 {T/cf 
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= V-e , ^ . (31) 
^/2(n + 2)(n + 8) 

As one lowers the temperature from the continuum high T limit, both R and 
Q vary as universal functions of r/T 1 '^ xv >: R decreases and Q increases as 
we lower the temperature into the magnetically ordered region (r < 0), while 
R increases and Q decreases as we lower the temperature into the quantum 
paramagnetic region (r > 0). 

The values in (|3l| ) are the key to the argument justifying the use of a 
classical dynamical model for small e. From (E8|) it is clear that the charac- 
teristic <& a fluctuations have an energy of order cy^R ~ \flT . As in Section |[ 
this energy is parametrically smaller than T, and so the occupation number 
of the relevant $ Q will be given by their classical equipartition value. This 
also means that the classical fluctuation dissipation theorem is obeyed, and 
the static susceptibility, x(k)> computed from ( p8|) by 

n 

T X (k) = -J2(\^(k)\ 2 ), (32) 

a=l 

is also the equal-time correlation of <j) a . 

We can now specify the recipe to compute the unequal time correlations 
of 4> a , in manner which parallels Section |3|. The fundamental Poisson bracket 
is 

{$a(aO, V-p{x')} PB = 5 a0 6{x - x'), (33) 

and the Hamiltonian H.2C then leads to the Hamilton-Jacobi equations of mo- 
tion 

^ = {® a (x),H 2C } PB 

= c 2 n Q , (34) 
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and 



— _ {n Q (x),H 2 c} PS 



at 

u 
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Vl$ a -R$ a - T (*|)#„, (35) 



The equations (f28|) , ( |34| ) and (|35|) define the central dynamical non-linear wave 
model of this section. We will compute correlations of the field <& Q at unequal 
times, averaged over the set of initial conditions specified by (^). Notice all 
the thermal 'noise' arises only in the random set of initial conditions. The 
subsequent time evolution is then completely deterministic, and precisely con- 
serves energy, momentum, and total O(n) charge. This should be contrasted 
with the classical-dynamical models studied in the theory of dynamic criti- 
cal phenomenali5li3, where there are explicit damping co-efficients, along with 
statistical noise terms, in the equations of motion. 

The dynamical model has been defined above in the continuum, and so we 
need to consider the nature of its short distance singularities. As in Section [| 
we assert M that the only short distance singularities are those already present 
in the equal time correlations analyzed earlier. These were removed by the 
simple renormalization in (p9|), which is therefore adequate also for the unequal 
time correlations. With this knowledge in hand, we can immediately write 
down the universal scaling form obeyed by the relaxation function by simple 



arguments based upon analysis of engineering dimensions. The analog of (25) 
is now 

1 _ / k LO 



KM = _*^_,_, e j, m 

where is a universal function we wish to determine. 

It now remains to solve the dynamical problem specified by ( p8|) , (^) and 
(|35|), and so determine ^n- For small e, the dimensionless strength of the non- 
linearity Q in ([^) is small; nevertheless we cannot use perturbation theory in 
Q, because this fails in the low frequency limit. In other words, the dynamical 
problem remains strongly coupled even for small e. 

The only remaining possibility is to numerically solve the strong-coupling 
dynamical problem. Formally, we are carrying out an e expansion, and so 
the numerical solution should be obtained for d just below 3. However, it 
is naturally much simpler to simulate directly in d = 2, which is also the 
dimensionality of physical interest. Therefore, the approach to the solution 
of the dynamic problem in the quantum critical region breaks down into two 
systematic steps: (i) Usa-the e — 3 — d expansion to derive an effective classical 
non-linear wave problem£3 characterized by the couplings R and Q . (ii) Obtain 
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Figure 6: Numerical resultsliiJ for thezero momentum scaling function ^75. (0, ui, Q) for the 
relaxation function, 1Z, appearing in (pa) for the d = 2, n = 1 model Zq. The dimensionless 
frequency ZU = lu/c^/R. Results are shown for Q = 25 (dots), Q = 30 (short dashes), 6 = 35 
(long dashes) and Q = 40 (full line). The high T limit value of Q in Q evaluates to Q = 35.5 
at e = 1 and n = 1. 



the exact numerical solutionis of the classical non- linear wave problem at these 
values of R and Q directly in d = 2. This division of the problem into two rather 
disjointed steps is also physically reasonable: it is primarily for the classical 
thermal fluctuations that the dimensionality d = 2 plays a special role, and 
the cases n — 1,2 (which have non-zero temperature phase transitions above 
the magnetically ordered phase) and n > 3 (which do not have a non-zero 
temperature phase transition) are strongly distinguished — so it is important 
to treat these exactly; on the other hand, the e = 3 — d expansion provides a 
reasonable treatment of the quantum fluctuations down to d = 2 for all n. 

Figs 6-8 contain the results of a recent numerical computation of the scal- 
ing function in (36) at k = 0. These results are the analog of Fig 4 for the 
fsing chain and Fig 5 for the d = 1, n = 3 case. They show a consistent trend 
from small values of Q and large values of n to large values of Q and small 
values of n, and we discuss the physical interpretation of the two limiting cases 
in turn. 

For smaller Q and larger n, we observe a peak in 7?.(0,w) at a non-zero 
frequency. This peak is the remnant of a delta function obtained in the large 
n limit at a frequency u> ~ T '. In the present computation, it is clear that 
the peak is due to amplitude fluctuations as $ Q oscillates about the mini- 
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Figure 7: As in Fig 6 but for d = 2, n = 2. The values of Q are now Q = 2Q_(short dashes), 
Q = 30 (long dashes) and Q = 40 (full line). The high T limit value of Q in (plh evaluates to 
g = 29.2 at e = 1 and n = 2. 




Figure 8: As in Fig 6 but for d = 2, n = 3. The values of Q are now Q = 2 (s hort dashes), 
Q = 25 (long dashes) and 5 = 30 (full line). The high T limit value of Q in (plh evaluates to 
g = 24.9 at £ = 1 and n = 3. 
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mum in its effective potential at $ Q = 0. As Q is reduced, we move out 
of the high T region into the low T region above the quantum paramagnet, 
and this finite frequency, amplitude fluctuation peak connects smoothly with 
the quantum paramagnetic quasiparticle peak. Of course, once we are in the 
quantum paramagnetic region, the wave oscillations get quantized, and the 
amplitude and width of the peak can no longer be computed by the present 
quasi-classical wave description — we need an approach which treats the excited 
particles quasi-classically. 

For larger Q and smaller n, the peak in 72.(0, to) shifts down to to = 0. 
The resulting spectrum is then closer to the exact solution for d = 1, n = 1 
presented in Fig 4. As Q increases further, the zero frequency peak becomes 
narrower and taller. How do we understand the dominance of this low fre- 
quency relaxation ? For n > 2, there is a natural direction for low energy 
motion of the order parameter: in angular or phase fluctuations of <1> Q in a 
region where the value of is non-zero. Of course, the fully renormal- 
ized effective potential controlling fluctuations of has a minimum only at 

= 0, as we are examining a region with no long range order. However, for 
these values of G, there is a significant intermediate length scale over which the 
local effective potential has a minimum at a ^ 0, and the predominant 
fluctuations of $ Q consist of a relaxational phase dynamics. 

The above reasoning has been for the cases with continuous symmetry, 
n > 2. However, closely related arguments can also be made for n = 1. In this 
case, in a region where |$ a | locally takes a non-zero value, there are low-energy 
modes corresponding to motions of domain walls between oppositely oriented 
magnetic phases. Indeed, precisely such a domain wall motion was mentioned 
for the d = 1, n = 1 case, and was argued to be behind the relaxational peak 
in Fig 4. 

Even in a region dominated by angular (or domain- wall) fluctuations about 
a locally non-zero value of |$ Q |, there could still be higher frequency ampli- 
tude fluctuations of \<f> a \ about its local potential minimum. This would be 
manifested by peaks in 7Z(Q, u) ) both at to = and at a non-zero frequency. 
A glance at Figs 6-8 shows that this never happens in a well-defined manner. 
However, for n = 1, we do observe a non-zero frequency shoulder in 1Z(0, lu) at 
Q = 35, along with a prominent peak at lu — 0: this indicates the simultaneous 
presence of domain wall relaxational dynamics and amplitude fluctuations in 
|$ Q |. Readers will also recognize the similarity of this with the shoulder in 
Fig 5 describing the high T limit of the d = 1, n = 3 case. For the other cases 
in d = 2, we do not see a clear signal of the concomitant amplitude and angular 
fluctuations: it appears, therefore, that once angular fluctuations appear with 
increasing Q, the non- linear couplings between the modes reduce the spectral 
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weight in the amplitude mode to a negligible amount. 

It is interesting to examine the above results at the value of high T limit 
for Q in ^ evaluated directly in e = 1. We find Q = 35.5, 29.2, 24.9 for 
n = 1, 2, 3, and these values are very close to the position where the crossover 
between the above behaviors occurs. The n = 1 case has a clear maximum 
in 1Z(0,lo) at u> = (along with a finite frequency shoulder), while there is a 
more clearly defined finite frequency peak for n — 3. 

In closing, we note that there is a passing resemblance between the above 
crossover in dynamical properties as a functimL pf £/, and a well-studied phe- 
nomenon in dissipative quantum mechanics E3li3 - EJ: the crossover from 'coher- 
ent oscillation' to 'incoherent relaxation' in a two-level system coupled to a 
heat bath . However, here we do not rely on an arbitrary heat bath of lin- 
ear oscillators, and the relaxational dynamics emerges on its own from the 
underlying Hamiltonian dynamics of an interacting many-body, quantum sys- 
tem. Our description of the crossover has been carried out in the context of a 
quasi-classical wave model here, but, as we noted earlier, the 'coherent' peak 
connects smoothly to the quasiparticle peak in low T paramagnetic region — 
here the wave oscillations get quantized into discrete lumps which must then 
be described by a 'dual' quasi-classical particle picture. 



5 Conclusions 

We have described the high temperature relaxational dynamics for a number 
of models in spatial dimensions d = 1,2. This dynamics is a property of a 
renormalizable, interacting continuum quantum field theory. Two cases can 
be further distinguished: 

(i) The excitations of the theory retain a non-zero scattering amplitude at high 
energies and temperatures: the models of Section || and^ are of this type. For 
these, the only characteristic energy scale controlling the density and interac- 
tion strength of the excitations becomes T itself, and so the phase coherence 
time, and the inverse relaxation rate, are universal numbers times Ti/ksT. As 
a result, quantum and thermal fluctuations contribute equally to the phase 
relaxation. (However, in Section ^ we did develop an expansion in which the 
universal prefactors of h/ksT became numerically large and so the long time 
relaxation was described by an effective classical model.) 

(ii) The theory becomes asymptotically free at high energies, and so the scat- 
tering amplitude of the excitations vanishes at large T. The model of Section]^ 
is of this class, and has a phase coherence time and inverse relaxation rate of 
order (ft/fe^T) ln(fcsT/ A), where A is an energy scale characterizing the low 
energy theory. These times are parametrically larger than Ti/kBT and so the 
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relaxational dynamics is classical. 

We have developed a fairly complete description of the dynamical corre- 
lations of these models, and our results should be testable in experiments on 
compounds the cuprate superconductor family, Heisenberg spin chains, and 
double layer quantum Hall systems EjO. 
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